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The broadband sound radiated by a subsonic rotor subjected to 
turbulence in the approach stream has been analyzed. The power 
spectral density of the sound intensity has been found to depend on a 
characteristic time scale — namely, the integral scale of the turbulence 
divided by the axial flow velocity — as well as several length-scale 
ratios. These consist of the ratio of the integral scale to the acoustic 
wavelength, rotor radius, and blade chord. Due to the simplified 
model chosen, only a limited number of cascade parameters appear. 
Limited comparisons with experimental data indicate good agree- 
ment with predicted values. 


The noise radiated by rotating machinery is known to consist of a 
broadband, random signal on which are superimposed a number of dis- 
crete peaks. These characteristics of the radiated noise spectrum have 
been observed on a wide range of turbomachinery such as aircraft com- 
pressors (refs. 1, 2, 3), turbines (ref. 4), fans (ref. 5), and marine pro- 
pulsors. 

The discrete peaks, which are due to mutual interference effects between 
blade rows, were of greatest concern as a cause of public annoyance and 
were the first to be investigated (refs. 6, 7, 8) . This source of sound is now 
well understood and its reduction has met with a certain measure of 
success (ref. 10). 

Broadband radiation has received relatively little attention in aircraft 
applications. The vibratory forces generated by turbulence, however, 
were of considerable interest in certain naval applications since marine 
propellers operate well within the boundary layer of a ship’s hull. In this 
case, the length scales of the turbulence and the propeller are of com- 
parable magnitudes, and, hence, strong thrust and torque fluctuations 
can arise (ref. 11). 

The random pressure fluctuations within a turbomachine which give 
rise to broadband radiation are due to the following phenomena: 
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(1) Turbulence in the approach stream generated, for instance, by 
upstream blade rows or by the combustion system. 

(2) Boundary layer turbulence at the duct walls. 

(3) Vortex shedding from the blade trailing edges. 

(4) Random variations of the tip vortex strength. 

Sharland (ref. 5), Smith and House (ref. 3), and Smith and Bushell (ref. 
4) conclude from their measurements that the major source of broadband, 
random noise in compressors and turbines is turbulence in the approach 
stream. This has also been found to be the case in marine propulsors 
(ref. 11), although vortex shedding from the blade trailing edges can — 
under certain conditions — become the dominant source over a range of 
frequencies. 

In this paper, the sound radiated from a rotor subjected to turbulence 
in the approach stream, at subsonic speeds, will be considered. 

CORRELATION OF PRESSURE FLUCTUATIONS 

Curie (ref. 12) and Doak (ref. 13) have shown that the dipole radiation 
from a turbulent flow containing rigid surfaces can lit' expressed in terms 
of the distribution of fluctuating pressures acting over the surfaces. In 
order to calculate the characteristics of these pressures, imagine the rotor 
blades to be subdivided into an arbitrarily large number of surface 
elements. The fluid velocity is measured with respect to a coordinate 
system a (a'= 1', 2', 3') that is fixed in space (fig. 1). An unprimed 
system a (a= 1, 2, 3) is fixed to the rotor and is so oriented that one of 
its axes coincides with the centerline of the rotor shaft. 

The time-dependent pressures acting on the various surface elements 
are related by virtue of spatial and temporal correlation of the velocity 
fluctuations in the approach stream as well as by virtue of induction 
effects that take place between adjacent blade elements. The acoustic 
radiation depends on the correlation area for the pressure-difference 
fluctuations; because of induction effects this area may differ considerably 
from the correlation area of the turbulence. 

Since the physical quantities involved in this problem are tensors, the 
index notation — including the summation convention in the case of 
repeated indices — constitutes the most convenient choice. Superscripts 
are used to denote the direction along one of the coordinate axes, whereas 
subscripts are used to denote the rotor blade element involved. For 
example, Uk B (r') denotes the component of the fluctuating velocity at time 
t in the direction /3 of the rotating reference frame at the blade element k. 
Similarly, F,*“ s (f,r') denotes the aerodynamic force acting on the fth 
blade element in the direction a at the instant of time t caused by a 
velocity fluctuation of unit magnitude in the direction 8 to which the Mh 
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Figure 1. — Coordinate system. 


blade element was subjected at the instant of time r'. Finally, (fit) 
indicates the aerodynamic force acting on the ith blade' element at time t 
in the direction a. In terms of these quantities, and neglecting higher order 
terms, the lift force is given by 


(fit) = 


f (<,/)«,*(/) dr' 

" QO 


( 1 ) 


where 


with 


a, 0 = 1 , 2 , 3 
i, k= 1, 2, . . . n 

r'<t 


In most cases, the aerodynamic force tensor is time-invariant, and 
equation (1) can be written as a convolution integral: 


■«)-/■ 


F ilc ^(t)u k ^t-r) dr 


( 2 ) 


where 


■ t t' > 0 


Since the aerodynamic forces (fit) are random functions, they are deter- 
mined statistically by the complete system of joint-probability dis- 
tributions of the values of the functions at any n values of t, where n may 
take any integral value. From an engineering standpoint, the correlation 
tensor ((i a (t)(f(t+r ) ) and its Fourier transform are the most significant 
quantities. In forming average values of the forces and fluid velocities, 
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we assume that the random processes are stationary and ergodic. The 
mathematical expectation 

1 t T 

Z?[t,“«K/(t+r)]= lim - / dt 

T- oo * J 0 [ 

= <(r) J (3) 


can be expressed in terms of the aerodynamic force functions and the 
velocity fluctuations in the form 


i$(r)= F*” (r,) dr, f FjfiHjt) dr* 

J a J a 

1 /■" 

X hm — / u k y (t — T )u, s (t+T — r 2 ) dt 

T-*o o * ■'0 

= { F ,t“' l '(r,) dr, f FjJ»(rt) dnR kr y> (r) 

n ^ n 


(4) 


Itkr yt (r) is the velocity correlation tensor for the points located at 
yt and y r , respectively. The spectrum tensor of the force fluctuations can 
now be found by taking the Fourier transform of the correlation tensor: 


(<*>) = - [ h,/ s (r)[exp tor] dr 

TT *'_oo 

= f T, t “' 1 '(ri)[exp 7«n] dr, f F [exp( — iW 2 )] dr -2 

* n ^ n 


(5) 


X 


J /?tr' |,s (r)Cexp( — nor)] dr 




where // (co) represents the aerodynamic frequency response function, 
Gkr yi {(F) is the spectrum tensor of the turbulence, and i = v 7- 1- 


APPLICATION TO A ROTOR 

As an illustration, consider a rotor placed in a flow containing homoge- 
neous and isotropic turbulence. In this case, the velocity correlation 
tensor can be expressed in terms of the distance r between two points in 
the flow field and the mean square value of the velocity fluctuations w 2 as 
follows : 
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Rkr aP (r) =u- 


rpOty%{i 

2 r 


3f 

dr 


+ /+ 


- r 
2 dr 


5 “0 


(6) 


where S a ^ is the Ivroenecker delta. The coefficient of longitudinal correla- 
tion/ can be approximated by an exponential function; namely 


/(r)= exp 



(7) 


where A/ is the integral scale of the turbulence. Furthermore, the velocity 
correlation tensor is approximated by 


Rkr ae (r)^ 



Ru^(c x t) 


( 8 ) 


where q = | y* — y r | is the distance between blade elements k and r. 

We now calculate the spectrum tensor of the force fluctuations normal 
to the plane of the rotor. In equation (5), a and /3 are both assigned the 
value 1, which corresponds to the direction normal to the rotor disk. The 
indices 7 and 5 take the values 1 and 2 only, since axis 3 coincides with 
the spanwise axis of a blade. From equation (8) the velocity correlation 
tensor has the components 




R kr u (r)=R k r 2 *(r)=0 

where c x is the axial velocity, as shown in figure 2. 


(9) 



u 



Figure 2. — Velocity diagram for 
an axial flow pump or com- 
pressor stage. 
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The Fourier transforms of these correlation functions can readily be 
calculated: 



The spectrum tensor of the force fluctuations therefore reduces to 

<?&(«) = [#«»(«) TlHjr" («) ](?*,” («) 

+ IHa 1 - («) ]*[//„“ («) ]<?*r 22 (co) (11) 


The simplest expression for the aerodynamic response function is that 
of Sears (ref. 14). This expression is based on two-dimensional, incom- 
pressible, thin-airfoil aerodynamics and ignores mutual interference 
between blade elements. The task of establishing aerodynamic response 
functions for cascades is yet to be accomplished. Recent experiments that 
have been performed at Pennsylvania State University, however, show 
that strip theory gives satisfactory predictions providing the value of the 
local lift slope is appropriately corrected. This correction consists in 
replacing the value of the lift slope — which is equal to 2ir in Sears’ equa- 
tion — by the steady-state local value. In terms of Sears function, the 
values of //(«) are given by 


where 


#,-,■“(«)=//<(«) sin 2 /3; 
7/,i 12 (co) = §//,( co) sin 2/3; 


( 12 ) 


//(co) 2-rrpwbSrK (k) 
K(k ) Sears function 

c ob 


K 


w 
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2b blade chord 

8r spanwise distance along the blade 

w relative fluid velocity 

/S angle between the flow relative to the blade and the axial 

direction, as shown in figure 2. 

We can further simplify our final expression if II (c o) is chosen as the 
mean value for the entire blade span. In this case 

=|Z/(«) \ 2 ^(2irw P by T ~-8r i 5r j 

ir ,'_1 y = i ( 1 I 2lT Kj 

(13) 

The spectrum of the total normal force driving the rotor can now be 
found from equation (11), using equations (10), (12), and (13). It has 
the form 


G(r) =*■ (pcfR 2 ) 



(14) 


where <t> = c x /U is the flow coefficient, B is the number of rotor blades, and 
u/c x is the turbulence level in the approach stream. 

The function S(R/A f ) represents a “correlation area” and is given by 




— (hi dr,-; 


r, 


(Id) 


The integration is performed from the hub radius r h to the tip radius 
R. This function has been plotted in figure 3 for rotors having various 
numbers of blades and hub-to-tip ratios. 

The frequency dependence of the total force driving the rotor is given 
by /( T) , which has the form 


/( r) = 



(10) 


where 


o)A / 

('X 


2b 

e=ir — cos 
A / 


0 


(7 = 2 
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Figure 3.— “Correlation area" in the spectrum function of the axial force. 


EXPERIMENTAL INVESTIGATIONS 

To verify the theory, an experiment was conducted in the water tunnel 
of the Ordnance Research Laboratory at Pennsylvania State University 
(ref. 11). This tunnel has a test section that is 4 ft in diameter and 14 ft in 
length. Velocities as high as 80 fps can be achieved, and the static pressure 
can be varied from 3 to 00 psia. The settling section of the tunnel contains 
a honeycomb of largo lcngth-to-diameter ratio that reduces the turbulence 
level in the test section to about 0.1 percent. 

A propeller was used for this investigation. The propeller had 10 blades 
with a constant chord length of 1 in. and a radius of 4 in. The design 
static thrust coefficient based on propeller disc area is 0.183, and the 
advance ratio at the design thrust coefficient is 1.17. The propeller and 
its installation in the water tunnel are shown in figure 4. 

A special balance was designed for measuring the unsteady thrust force 
of the propeller. The arrangement used is shown in figure 5. A piezo- 
electric crystal is mounted in a steel cup at the end of the propeller shaft. 
After assembly, the cup is positioned by means of set screws until the 
hemispherical ball bundl'd to the crystal lies on the exact center line of the 
shaft, thus minimizing the crystal’s response to bending distortions of the 
shaft caused by hydrodynamic moments acting on the propeller. The 
frequency response and the linearity of the balance are shown in figure 6. 

During the tests, turbulence was generated by means of grids mounted 
20 mesh sizes upstream of the propeller. Two grids were used with mesh 
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Figure 6. — Frequency response and linearity of balance. 


sizes of 4 in. and 6 in., respectively. The first had a solidity of 0.34 and 
was fabricated with f-in. -diameter rods. The second had a solidity of 
0.27 and consisted of |-in. -diameter rods. Figure 7 shows a view of the 
propeller taken from a position upstream of the 4-in. grid in the water 
tunnel. 

The measurements are compared in figure 8 with theoretical predictions. 
The “humps” in the measured spectrum occur at blade passage frequency. 
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Figure 7. — Downstream view of propeller through the J,-inch grid. 


SOUND RADIATION FROM A ROTOR 

The statistical measure of greatest engineering interest is the power 
spectral density of the sound pressure fluctuations. This is related to the 
power spectral density of the force fluctuations by a Green’s function 
g(x, y;a>) as follows: 

P(x;w)= (17) 

*=i ;= i 

If the rotor is housed in a duct, the radiation into the free field will be 
affected by the characteristics of sound propagation in the duct and by 
its intake geometry. These effects have been considered by Tyler and 
Sofrin (ref. 6), as well as by Morfey (ref. 15), for pure tone radiation. 

In order to simplify the present problem, the Green’s function chosen 
will be that of a dipole whose axis is normal to the rotor plane. Comparison 
of predicted acoustic power with experimental measurements on ducted 
fans will still be possible in cases where the dimension of the duct is com- 
parable with, or greater than, the quarter wavelength of the sound 
generated. The appropriate Green’s function for the far field is given by 

xtc cos 0 

g(x,yj]u) = {expi'/cr} {expl-i/cr, cos (<?,-<?) sinfl]} (18) 

4irr 



dB re 1 lb in a 1-Hz bandwidth dB re 1 ,b in a 1_Hz bandwidth 
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TEST No. 5531, RUN No. 2 

DISTANCE BETWEEN GRID AND PROPELLER - 80 INCHES 
WATER-TUNNEL TURBULENCE LEVELWITHOUT GRID ■ 0.0011 U 
TURBULENCE LEVEL AT PROPELLER DUE TO THE GRID - 0.03 U 
TUNNEL VELOCITY ■ 15.4 ft/sec 
PROPELLER ADVANCE RATIO - 1.22 



0 100 200 300 400 500 600 


FREQUENCY IHzl 
a. — ^-inch grid mesh. 


TEST No. 5530, RUN No. 2 

DISTANCE BETWEEN GRID AND PROPELLER - 120 INCHES 
WATER-TUNNEL TURBULENCE LEVEL WITHOUT GRID - 0.001 U 
TURBULENCE LEVEL AT PROPELLER DUE TO THE GRID - 0 03 U 
TUNNEL VELOCITY - 15. 1 ft/sec 
PROPELLER ADVANCE RATIO - 1. 22 



FREQUENCY (Hz I 


1). — 6-inch grid mesh. 


Figurk 8. — Power spectral density of the. propeller response. 
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where k = w/c. The polar coordinates used in this equation are defined 
in figure 9. 

A further restriction is imposed by the requirement that the eddies be 
compact on a wavelength scale so that retarded time changes will be 
negligible. This implies that the acoustic wave length X must be much 
larger than the integral scale A/. 

The spectral density of the sound intensity at x follows from equations 
(17) and (18) : 



(19) 

where M = c x /c is the axial flow Mach number and g(T) = F/(r). 

In deriving the function F, it is assumed that the solidity of the rotor 
is such that it can be treated acoustically as a disk. This implies that the 
acoustic wavelength must be much larger than the blade chord and the 
blade spacing. Under these conditions F has the form 



X [ Jo{kRQ sin 6) -\-Jo(kRZ sin 6) ]r,- dr,- d<pi rj drj d<pj 

(20) 

where J 0 denotes Bessel functions and r/, denotes the hub radius. The 
quantities Q and Z are measured on the rotor disk and are given by 

Q = ^ (r< cos tpi—Tj cos <pj) 



z 


Figure 9. — Polar coordinate system. 
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z = 


1 

R 


(ri sin ipi—Tj sin $,) 


Finally, the intensity can be integrated over a large spherical surface 
to yield the spectral density of the sound power : 


dW 

dr 




( 21 ) 


where 


F ‘(j/ kR )’ C F ( e -t/ kR ) a 


sin 6 cos 2 0 dd 


The function F\ has been calculated numerically and is plotted in 
figure 10 for a range of values of kR and R/ A/. 

As an example, the theory has been applied to the following case: 


Rotor hub-to-tip radius ratio 


n 

R 


= 0.5 



Figure 10 .—Radiation function in the spectral density of the sound power. 
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Flow coefficient 
Axial flow Mach number 


<£ = 0.64 
M = 0.25 

Ratio of integral scale to tip radius — = 0.5 

R 

This choice of A/ corresponds to the span of the blades. The resulting 
spectrum is shown in figure 11. The vertical scale is proportional to 
T(dW/dT); i.e., a constant-percentage bandwidth, and the horizontal 
scale consists of the ratio of the integral scale A/ to the sound wavelength 
X. The peak of the spectrum occurs at a value of A//X = 0.16. The pre- 
dicted spectrum level becomes inaccurate at values of A//X>0.25 since 
the integral scale of the turbulence and the wavelength of the sound 
become of comparable magnitude. A much more rapid, exponential 
decrease in the value of Y(dW/dT) than indicated in figure 11 would 
be expected. 

The directivity of the sound pressure is illustrated in figure 12, which 
shows that the sound becomes progressively less directional as the integral 
scale A/ decreases. 

Af/X 

0 0.008 0016 0032 0063 0.125 0250 0500 

rd# 
dr 

-10 

db 

-20 

Figure 11 . — Broadband sound-power spectrum. 
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R/Af = 0 

R/Af = 4 0 


Figure 12 . — Directivity of sound 
radiation. 
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SUMMARY AND CONCLUSIONS 

Broadband sound radiation from turbomachinery appears to be mainly 
due to turbulence in the approach stream. This turbulence may be gen- 
erated by one of several mechanisms such as upstream blade rows or, 
in the case of turbines, the combustion system. 

As a preliminary step in the calculation of the power spectral density of 
the radiated sound pressure, the spectrum tensor of the axial force fluc- 
tuations acting on a rotor have been found. These have been expressed in 
terms of aerodynamic response functions and the correlation tensor of the 
turbulence. At the present time, there is little information concerning the 
characteristics of the turbulence in turbomachinery. Also, expressions for 
the aerodynamic response functions are limited to two-dimensional, 
incompressible, thin-airfoil theory. Research in both areas is desirable. 
In order to make further progress and establish the general characteristics 
of the radiated sound, a simple case was considered. The turbulence in 
the approach stream was chosen to be isotropic and homogeneous, and 
the spectrum of the total force' normal to a rotor disk was calculated. The 
predicted spectrum levels were compared with experimental results 
obtained on a 10-bladed propeller subjected to grid turbulence. Good 
agreement was observed. 

The power spectral density of the sound intensity and of the sound 
power were next calculated. The spectrum depends on a number of 
parameters, such as the turbulence level, a characteristic time scale con- 
sisting of the ratio of the integral scale of the turbulence divided by the 
axial flow velocity, and characteristic length scales such as the ratio of 
the integral scale of the turbulence to the acoustic wavelength, the radius 
of the rotor, and the blade chord. 

Due to the simplicity of the model used, only a limited number of 
cascade parameters enter the equations. In a more complete description, 
additional cascade' parameters — such as the solidity and stagger angle — 
would also appear. As an example, the broadband sound-power spectrum 
for a turbomachinery rotor was calculated and plotted in figure 11. A 
flow coefficient of 0.04, a hub-to-tip ratio of 0.50, and an axial flow Mach 
number of 0.25 were chosen. The broadband sound power radiated has 
been estimated at 140 dB re 10~ 12 watts for a turbulence level of 5 per- 
cent of the axial velocity. 

LIST OF SYMBOLS 

B number of blades on rotor 

2 b blade chord 

C cascade parameter defined in equation (16) 
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c speed of sound 

c x axial flow velocity 

e cascade parameter defined in equation (16) 

F aerodynamic force function, defined in equation (1) 

/ coefficient of longitudinal correlation of the turbulence 

G(u) spectrum function 

H (w) aerodynamic frequency response function 

K(k) Sears function 

k acoustic wave number 

( local lift force 

M Mach number of axial flow 

n number of blade elements 

Q,Z projections of q on the x and y axes, respectively 
q distance between blade elements 

R tip radius of rotor 

r distance, as variously defined 

r h hub radius of rotor 

t,T time 

V tangential velocity of rotor 

u fluid velocity 

W sound power 

w flow velocity relative to the blades 

x,y position vectors 

Greek Letters 

(i angle between the flow relative' to the blade and the axial direction 

T frequency parameter, defined in equation (16) 

$ correlation function 

< p flow coefficient, defined in equation (14) 

k reduced frequency parameter in Sears function 

A / integral length scale of the turbulence 

X wavelength of the acoustic radiation 

r time 

6,<p position angles in polar coefficient 

u > angular frequency 



APPENDIX: TABLE OF RADIATION FUNCTION 1 
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This table applies to a hub/tip ratio of 0.50. 
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DISCUSSION 


R. MANI (University of Massachusetts): It would seem from your 
description of the experimental apparatus that you have turbulence scales 
which are quite large compared to the propulsor size. It has been shown 
that you should get a predominance of blade-passage tones. If the eddies 
are really of substantial size, I do not really see how you would only get a 
broadband spectrum. 

SEVIK (author) : I believe that you are referring to a paper by Ffowcs 
Williams and Hawkings in which they show that blade-passage tones will 
also occur when the inflow is turbulent. They indicate that the bandwidth 
of the spectrum will be inversely proportional to the number of blades 
which are subject to the distortion field of a single eddy. Our test results 
do indeed show a “bump” around blade-passage frequency, which corre- 
sponds to their prediction. The bandwidth of the bump is about 40 cps, 
in agreement with their expression. Unfortunately, we also discovered 
after the tests that we had a blade resonance in the same frequency range 
and hence that there were two possible contributors to the energy. I 
agree that my theoretical approach does not predict the tones, but only 
the broadband energy. 



